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COMPENSATING POLARIZATION MODE DISPERSION IN FIBER OPTIC 

TRANSMISSION SYSTEMS 

RELATED APPLICATIONS 

This application is a continuation-in-part of U.S. Patent Application Serial No. 
5 09/505,313, filed February 16, 2000 and under 35 USC §1 19(e)(1), claims the benefit of U.S. 
Provisional Application Serial No. 60/243,522, filed October 23, 2000, which is incorporated 
herein by reference. 

BACKGROUND OF THE INVENTION 

This invention relates generally to fiber optics and more specifically to an apparatus 
10 for reducing the polarization mode dispersion in a fiber optic transmission system. 

Polarization mode dispersion (PMD) refers to variations in the time delay of a polarized 
optical signal traveling through an optical transmission system, such as a single-mode optical 
fiber. PMD arises in an optical fiber as a result of asymmetries in the optical fiber's core, 
such as core ellipticity created during optical manufacturing and bending stresses formed 
15 during installation. The asymmetries of the fiber's core cause random changes in the state of 
polarization (SOP) of optical signals propagating through the fiber. Different SOPs 
propagate through the optical fiber core at different relative speeds, e.g., some SOPs ti-avel 
faster and some travel slower, resulting in a pulse width distortion of a transmitted optical 
signal relative to an input optical signal. Additionally, the asymmetiies of the fiber's core are 
20 highly susceptible to enviromnental fluctuations, such as temperature or movement of the 
fiber, which occur as fast as milliseconds and resuh in a time varying pulse widtii distortion 

of the transmitted optical signal. 

The varying pulse width distortion is mathematically represented by a time delay 
between two orthogonally polarized principal states of polarization ("PSP") which form a 

25 convenient basis set to describe and characterize each SOP, and evaluate the effects of PMD 
in the fiber. Using tiie PSPs as a basis set, each SOP propagating tiirough an optical fiber is 
represented by a linear combination of the two orthogonally polarized PSPs. The varying 
pulse width distortion of the SOP is a fimction of a varying delay between the PSPs. 
Theoretically, each PSP experiences a time of flight difference through the optical fiber, 

30 commonly known as differential group delay, resulting in a time delay between tiie two PSPs 
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at fte fiber output. The output SOP is represented by a linear combination of the PSPs which 
are time delayed with respect to each other. A greater time delay between the PSPs 
corresponds to a larger relative difference between the input SOP pulse width and the output 
SOP pulse width. See tor example CD. Poole and R.E. Wagner, "A Phenomenologrcal 
Approach to Polarization Dispersion in lx,ng Single-Mode Fib^." Elect^nic Ut^rs. Vol. 
22 pp 1029-1030, September 1986, which is mcorporated by reference herem. 

' Optical fibers have a differential group delay (DGD) between the two PSPs on dte 
order of 0.1 psfcn. In older fiber optic cables, such as the cables used in terrestrial ne^vorks, 
the DGD is on the order of 2.0 psfcn and results in time delays of about 50 picoseconds for 
.remission distances of oriy several hundred kilometers. As the demand for faster ophcal 
data transmission increases, such as from gigabits per second to terabits per second, opttcal 
pulse width distorUon due to PMD will become one of the factors limiting data transmrsston 
rate. 

SUMMARY OF THE INVENTION 

A real-time optical corrxpensating apparatus reduces first-order PMD in an optical 
fiber by determining the PSPsofthe optical fiber and delaying one PSPwithrespect to the 

other. . 1 J- 

in one aspect, the invention features an optical compensating apparatirs melud>ng a 

polarization controller configured to receive an optical signal propagating through an optical 

, medium, determine the principal states of polarization of tire optical medium, and transform 

the optical properties of ti^e optical signal based on flte principd states of polarization. 

Embodiments of tins aspect may include one or more of tire followmg. The 

polarization controller includes a polarimeter configured to teceive and to measure the state 

of polarization of tite optical signal, =u,d a polarization transfonner arranged in the patit o th 

5 optical signal after the polarimeter and configured to tnnsfomt the polarization of the opttcal 

signd The polarization controller includes a central processing unit programmed to momtor 

the time averaged state of polarization of the optica signal and determine tite principle «ates 

of polarizationof tite optical medium based on the time averaged state of polanzatton. The 

cential processing unit analyzes time averaged sm.es of polarization via ttte relationsh.p. 
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_ , to determine the principle states of polarization by calculating 



the vector a . The central processing unit is programmed to determine a magnitude of time 
delay between the principle states of polarization based on the time averaged state of 
polarization. The central processing unit is programmed to determine the relative order, m 
5 time of the principle states of polarization. The compensating apparatus includes a delay 
controller arranged in an optical path of the optical signal after the polarization controller. 
The delay controller is configured to receive the transformed optical signal and to generate a 
signal proportional to the PMD time delay between the principle states of polarization. Tlie 
delay controller includes a DOP sensor, such as a polarimeter or spectrum analyzer. The 
1 0 delay controller includes a delay transformer arranged in the path of the transformed optical 
signal before the polarimeter. The compensating apparatus further includes a PMD 
generating module arranged in the optical path of the medium before the polarization 
controller and configured to receive the optical signal from the optical medium and generate 
a known PMD. THe compensating apparatus fiirther includes a calibrating unit arranged m 
15 anopticalpathoftheopticalsignalafterthepolarizationcontroller. The calibrating unit is 
configured to receive the transformed optical signal and determine the transformation of the 
polarization controller. The calibrating unit includes at least one detector for detecting one or 
more components of the optical signal. 

In another aspect, the invention features, an optical compensating apparatus for 
20 reducing PMD in an optical signal transmitted through an optical medium. The optical 
compensating apparatus includes a polarization controller configured to receive an optical 
signal propagating through the optical medium, determine the principal states of polarization 
of the optical medium and a magnitude of time delay between the principal states of 
polarization; and a delay controller arranged in an optical path of the medium after the 
25 polarization controller and configured to receive the transformed optical signal, wherem the 
polarization controller transforms the optical signal based on the principal states of 
polarization and the delay controller compensates PMD of the optical signal by reducing the 
time delay between the principle states of polarization based on the magnitude of time delay. 
Embodiments of this aspect may include one or more of the following. The 
30 polarization controller includes a polarimeter. The polarization controller includes a 
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polarization .ransfonne. ammged in path of the optical .ignal after the polarimeter. The 
delay controller is configured to generate a signal proportional to the PMD time delay 
between the principle states of polarization. TTte delay cont„>ller includes a DOP sensor, 
such as a polarimeter or spectrum analyzer. The delay controller includes a delay 
, transformer arranged in fte path of the transformed optical signal before the polarimeter of 
.he delay co„»ller. The delay controller includes a delay transformer arranged m the path 
ofthe transformed optical signal after the polarimeterofthe polarization controller and 

before a polarimeter of the delay controller. 

In another aspect, the invention features an optical compensating appamtus for 
,0 ^ducingPMDinanopticalsignaltransmittedthroughanopticalmedium. Theapparatus 
includes a polarization module configure to receive an optical signal propagating through 
d,e optical medium, determine the principal states of polarization ofthe optical meditun. and 
generate a s.gnal for transforming the polarization ofthe optical signal; and an opucal 
transformer arranged in an optical path of the medium after the polarization module arrd 
,5 configured to transform the optical signal and .educe a Ume delay between the pnncple 
states of polarization based on the signal received ftom the polarization module. 

Embodiments of this aspect may include one or more ofthe foUowmg. The 
polarization module includes a polarimeter. m optica, transformer includes a polarrzatton 
transformer and a delay controller. The polarization transformer and a delay controller each 
ao include a polarimeter. The polarimeter ofthe delay controller is configured to generate a 

signd proportional to the PMD time delay. The delay controller includes a delay transformer 
arranged in the path ofthe transformed optical signal before the polarimeter of *e delay 
controller and after the polarization transformer. 

,n another aspect, the invention feam^s a method of r^ucing PMD of an ophcal 
signalpropagatinginanopticalmedium. The method includes determining the princtpal 
state of polarizations ofthe optical medium with a polarization controller, and transfonmng 
the polarization of Ute optical si^ with a polarization transforming device based on the 
polarization of the first principal state of polarization. 

Embodiments of this aspect may include one or more of the following. The method 
30 further includes determining the time delay between a fnst principal state of polarizafion and 
a second principal state of polarization. The method fitrther includes delaying the first 
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Principe .a. of polaHzation wia. respec, .0 a second principal s.. of polarizaUo. 
LJdf«heri„c.ud=sdc.ayi«gaflrs,pri„dpals.,eofpoIari.atio„w,U>re.pec tea 

principal .a. of poiarization. T.e ™=d,od furihcr includes d=tc« U-e .la .ve 
:"n^i»e^f.eprinciplcs..csofpolari..on.Tl.e.ed,„d.^.«^^^^ 

a time averaged state of polarization of an optical signal ptopagatmg through the opt cal 
1 L. Tl^^stepof deterntiningtheprinciplestatesof polarizationincludesa^^^^^^ 
.i„,eaverageds.ateofpo,=uization. Thetinte averaged statesofpolartzattonfollowthe 

1 ? 1- , and determining the principle states of 

relationship, DOP = 1 '^^ I " ] _ (5 • f , )a 1 

. . .u The vector a is determined by analyzing rk 

nolarization includes calculatmg the vector a . The vector , . , , 

:lo..earcurvef.tUng.echni,ucsorn.atri.n.e.hods.T.e.nethod«^.^^^^^^^ 

calibrating a polarization controller for transforming the optical stgnal to d_^ 
plurality of control settings for the polarization controller by monttonng whteh control 
gTals when applied to the polarization cont.1 transform an input optical st^al havtng a 
p::cstate:;iarizaUonintoanou.putsi^.hav.ngditfer.tstateofpo^^^^^^^^ 

,„ another aspect, the invention features a ntettod of compensafng PMD of a„ opUcal 

signal propagating through an optical medium. The method includes recordmg the tunc 
Ig d Jeofpolarizationofanopticalsignalpropagating^roughanopttcalm^mm 

an. lyzing the time averaged state of ^larization of the optical st^a, vta a pertubattve 

cpansionof the time averaged state of polarization. 7 . to extract a plurality of parameters 
, rL.imeaveragedstateofpolarization.hatcharacterizethePMDof.heopt,cals.^, 

wherein the pertubative expansion is expressed as , . f„, 

DOP^ .\rMs„f H.-A<»^.r;.»)l?.J=-0l--')- The plurality of parameters, for 

example can be used to cha«cterize the location otthe principle states of polarization on me 
:™^s;el,themagnitude„ftimede,aybetwee„thepnneiplesta.esofpo— ,and 



25 



*^"':i:::r:— esinclud.ng,bu.notlimi.edto,oneormoreof.^^ 
followir.TheapparatusforcompensatingPMDoperate.inrea,rimeanddoesnotre,u,rea 



fast detector. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 isablock diagram ofaPMD compensating apparatus; 

Fig. 2 is a block diagram of a polarimeter shown in Fig. 1; 

Fig. 3 is a block diagram of a delay module shown in Fig. 1; 

Fig. 4 is a schematic representation of a Poincare sphere; 

Fig. 5A is a schematic representation of an optical signal propagating through an 

optical fiber of Fig. 1; 

Fig. 5B is a graphical representation of the optical signal propagating through an 

optical fiber of Fig. 1; 

Fig 6 is a schematic flow chart of the PSP algorithm; 

Fig. 7A is a schematic representation of the principal states of polarization of an 
optical signal propagating through an optical fiber of Fig. 1; 

Fig 7B is a graphical representation of the degree of polarization as a fianction of the 
time delay, x, between and the relative power, a, of each principal state of polarization; 

Fig. 8 is a cross-sectional view of the Poincare sphere of Fig. 4 taken along the 
equator; 

Fig. 9 is a block diagram of another PMD compensating apparatus; 

Fig 10 is a block diagram of another PMD compensating apparatus; 

Fig. 11 isablock diagramofaapparatus for calibratingapolarization transformer; 

Fig. 12 is block diagram of an expanded view of a delay module integrated with a 

calibration unit. 

DESCRIPTION OF THE EMBODIMENTS 

Referring to Fig. 1, compensating apparatus 10 includes a polarization controller 100 
and a delay controller 200. Compensating apparatus 10. when placed between an output 21 
of an optical fiber 22 and an input 235 of optical receiver 240. reduces the PMD of optrcal 
signals transmitted by optical transmitter 15 though optical fiber 22. 

Polarization controller 100 includes a lens 104. abeam splitter 105. a polarimeter 
no andapoiarizationt^nsformerlOS. Uns 104 positioned a. an input 102 of polarizatron 
comroller 100 coUimates optical signals (not shown) from output 21 of optical fiber 22 along 
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an optical path 160. Optical path 160 extend, ftotn inpt« end 102, through beam spliner 
105 polarization controller 108. and out output end 103. Beam splitter 105 redirects a 
portion of the optical signal propagating along beam path 160 into polarimeter 1 10 wh,ch 
detects the redirected optical signals and sends a series of electronic signals through cables 
122 to an I/O port 1 19 of a computer 120. Computer 120 uses the electronic signals m an 
algorithm stored in Ure computer's CPU to determine the principal states of polartzatton 
(PSPs) of optical fiber 22 and sends control signals to modify the setdngs of a first retarder 
140 and a second retaider 150 in polarization transfomter 108. First retarder 140, e.g., a 
quarter-waveplate, and second retarder 1 50, e.g., a half-waveplate, transform the polanzatton 
of the PSP such that light exiting polarization controller 100 is linearly polarized and altgned 
to the X- and y- axis of delay controller 20O. 

Delay controller 200 includes a delay module 170, a beam splitter 165, a polartmeter 
210 a controller 220, and amirror 202. An optical beam path 162 extends between an tnput 
16l' atrough delay module 170, and beam splitter 165. Mitror 202 reflects optical signals 
out of delay controller 200 through output 163 and into input 235 of receiver 240. Beam path 
162 a, input 161 is colUnear with beam path 160 from polarization controller 100 such dtat 
coUimated optical signals exiting output end 103 propagate along beam path 162. 

After the optical signals pass through delay module 170, beam splitter 165 redtrects a 
portion of the optical signal propagating along beam path 162 into polarimeter 210. 
Polarimeter 210 detects the redirected optical signals and sends a series of electronic signals 
via cables 21 5 to a control cUcui, 220. Polarimeter 210 is similar in s,ruc«« to polanmeter 
110 described below. Control circuit 220 uses flte electronic signals sent from polanmeter 
2,0 to determine the time delay between the PSPs in optical fiber 22 and then sends a control 
sigral via cable 1 30 to delay controller 170. The control signal modifies tite settings of delay 
controller 170 such flra, the time delay is reduced between the two PSPs transmitted duough 

outlet 2 1 of optical fiber 22. . , , , „ j 

Referring to Fig. 2, polarimeter 110 includes three beam splitters 1 14, 1 16, 1 17, and a 
mirror 119spaced along an optica beampath 112. Beamsplitters 114, 116, IH.and mirror 
1 1 9 couple optical signals propagating along beam path 1 12 u,™rds detector modules 1 14a, 
, U6a 117a, 119a, respectively. Detector module 114a includes a detector 114b for 

nreasuring the total power of an optical signal. Detector module 1 16a includes a polartztng 
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beam splitter U 6b and a detector assembly U6c having a first detector U6d and a second 
detector 1 16e. Similarly, detector module 1 17a includes a polarizing beam splitter llTb and 
a detector assembly 1 17c. Detector module 1 19a includes a polarizer 1 19f. e.g., a quarter- 
waveplate, a polarizing beam splitter 1 19b, and a detector assembly 1 19c. Each detector 
module measures specific optical properties of the optical signal and sends an electromc 
signal proportional to each measured property to computer 120 via cables 122. 

Referring to Fig. 3, delay module 170 includes an input polarizing beam splitter 171, 
an optical delay assembly 174, and an output polarizer 172. Polarizing beam splitter 171 
separates the two PSPs transmittal through optical fiber 22 and polarization controller 100 
such that PSPl, delayed with reject to PSP2, propagates along a fixed optical path 173 to 
polarizing beam splitter 172, and PSP2 propagates along a variable optical path 175. 
Variable optical path 175 includes optical delay assembly 174. e.g., a tianslatable (Arrows) 
mirror, which delays PSP2 with respect to PSPl. PSPl and PSP2 r^mbine in polar,z,ng 
beam splitter 1 72 and continue propagating along beam path 162. 

Referring to Fig. 4, a convenient and intuitive graphical representation of SOPs is a 
Poincar^ sphere 500. A SOP is defined in terms of a Stokes vector of Formula 1: 



20 
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30 



{^0 s, s, s,) 



(1) 



where 



51 =E\ - E y 

52 = 2ExEyCOs(0) 

53 = 2ExEySin(0) 

and B. and E, ar^ the magnitiKles of the x and y component electiic field complex 



(2) 
(3) 
(4) 
(5) 



amplitudes, respectively. 0 is the relative phase between Ute two. The parameters s„ s. and 
s, are represented by the relationship s.= (S,)/(S.), where i is 1, 2, or 3, and can be used to 
convert ti,e Stokes parameters to corresponding x, y. and z components m a three 
dimensional cartesian coordinate system. ,„ a fl»ee dimensional Cartesian coord,uate 
system, Poinoar6 sphere 500 is defined by a set of points containing all posstble SOPs. As 
seen in Fig. 4, all Uneax polarization smtes are located on an equator 5 10 of sphere 500, wtale 
left and right circular polarizations are located at a north pole 520 and a south pole 530, 

— 3— 



ni 



rl 20 



^^ocketNo.: 01997-296001 

respecUvely. All ote points represent elliptical polarizadons whieh lie away from equator 
510 and poles 520, 530. Each SOP on sphere 500 is identifiable by its latin.de 2w and 

longitude Tk by using equations: 

si = cos(2w)cos(2A,) 

(7) 

5 S2 = cos(2w)sin(2A-) 

S3 = sin(2w) 

where any two orthogonal SOPs lie directly opposite each other, e.g., a linear vertical 
polarization 512 is 180 degrees away from a linear horizontal polarization 514 on equator 

510. . J , , 

In operation, transmitter 15 sends a polarized optical signal to polarization modulator 

20 which modulates the state of polarization ("SOP") of the optical signal, e.g., from vertical 

to right circular to elliptical, with a frequency of about 10 kHz to about 100 MHz. The 

modulating frequency is fast enough to measure and compensate varying PMD on a 

millisecond timescale. The modulating frequency is limited by the response time of the 

16 detectors used in the polarimeters. 

Referring to Fig. 5A, as an input SOP 550 propagates from polarization modulator 
(not shown) through optical fiber 22 towards output21, the SOP of the signal randomly 
changes. Each SOP propagates through the optical fiber at a different speed, e.g., some SOPs 
travel faster and some travel slower, resulting in a varying pulse widfl> distortion of an opttcal 
signal 560 at output 21. To a first order approximation, optical fiber 22 has two discrete 
group delays, one for each of two orttrogonal PSPs, i.e.. PSP 562 and PSP 564. Referring to 
Fig 5B optical signd560,aSOP.isalinearcombinationofPSP 562 and PSP 564. Atime 
delay 565, e.g., 40 ps, between PSP 562 and PSP 564 creates PMD in the output signal, i.e., 
the width of signal 560 is greater than signal 550. 
25 At outtet 21, optical signal 560 propagates into compensating apparatus 10 which 

reduces tire PMD in signal 560 witi. polarization controller 100 and delay controller 200. 
Polarization contioller determines the PSPs of flte optical fiber 22 and tiansforms tire two 
PSPs to X and y linearly polarized states aligned with the x and y optical axis of delay module 
170. Delay contioller 200 measures tire time difference between the two tiansformed PSPs 
30 and delays one PSP relative to tire Other. 
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After exWng optical fibe, 22, optical signal 560 travels through polarization 
controller .00 along beam path 160, mttil bean, splitter 105 r^irects a portion, e g., about 
1% of the optical signal into polarimeter 110 for ^ysis. Tlte amount of opncal stgnal 
redirected into polarimeter is sufficient such fltat the redirected opUeal sig^l ts measurable 
by each of the detectors in polarimeter 110. 

ReferringbaclttoRg.2,beamsp„t,ers 1.4, ,16, 117 everiy divtdetheopUca. stgnal 

entering po.arimeter 1 .0 into four separate opticM signals. Detector module 1 14a measures 
,he power ofthe firs, optical si^U.e..S.. Detector module .16a analyses the second 
signal by measuring the difference between opUcal signals having polarization componen^ 
„ entedinthexandydirection,i.e.,S,=EVEV optical signalshavingo^y...^^^ 

polarization components result in a measurement of .1 by detector module 1 16a, ^ opttca. 
signals havingonly y oriented polarization components result in a measurement of -1. 
4ctormod„lell7aana.yzesthethirds.^bymeasuring.hedifferencebetweenop.,cal 

sign^s having polarization components oriented 45 degrees with respect « 

direction, i.e., S. - E^.« = 2E.E^s(0) (where 0 is the phase between E, and E,) Dete tor 

module 1 1 9a analyzes the fourth signal by measuring the difference betwe^ opt^ s^als 

having left and right circular polarization components, i.e.. S, - 2E,E,sin (0) where 0 .s the 

phase between E. and E,. Optical signals having only right circular polarization results m a 

measurements of+1. and left circular polarization results in a measurement of -1^ 

;rmeterU0measureseachof.heS,oltesparame.ers.E,uations2-5,andsendselectr„ntc 

signals proportional to each measurement to computer 120. 

Lfling.oFig.6.compu.er.20receives,heelec.ricalsignalsftompolanme«m 

andrunsanalgorithm 600 stored inthecomputer'sCPUtodetermine^e^ 

on aPoincar* sphere (S.O), to calculate a de^ of polarization (DOP, for ea^h SOP (S20). 

. .0 determine the two orthogonal PSPs for the opUca, flber (S30). and to ^^^^^^ 

polarization transformation which converts the PSPs of the fiber to hnear x and y polartzatton 
states aligned with the x and y optical axis of delay module 170 (S40). 

The CPU determines the location of the SOP on the Poincare sphete by relating each 
of tite electiical signals ftom tite polarimeter to its corresponding Stokes parameter, equations 

,0 6-8 and then calculating w and X. i.e.. tite coordinates of the SOP on the Poincare sphere. 
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Each of the Stokes parameters is also used by the CPU to calculate degree of polarization 
(DOP). The ratio of Stokes parameters shown below 



DOP" = i 



(9) 



is less than or equal to unity. 

Referring to Fig. 7A, two PSPs, i.e.. square shaped pulses of leugth ... have a first 
order PN© delay of x and a relative power with respect to each other described by the 
parameter a. To a first order approximation, i.e., when x/x. is less than 1, the DOP of the 
two square PSPs, equation 9, can be rewri«en as a ftnCion of delay, t. and relaUve power, a: 



1/2 

DOP(T,a) - [1 + 4 a (t/xo) ((x'xo)-2)(l- a)] ■ 



(10) 



,0 A plot of Equation 10 (Fig. 7B) graphic^ly shows how DOP depends on delay and relattve 
power of the PSPs. At constant delay, x, DOP is at a minimum when both PSPs have equal 
power whereas DOP is at a maximum when only one PSP has all of the power, i.e., a ,s 1 or 
0 respectively. At constant relative power, a, DOP is inversely related to the delay between 
fte two PSPs. A, a constant time delay between the two PSPs, the DOP depends on a>e SOP 
,5 When fl.e SOP is a 50/50 mixtirre of both PSPs, i.e., each PSP has equal power, the DOP wtU 
be at a minimum, whereas the DOP will be unity when the SOP is aligned witi, one PSP, >.e., 
one PSP has all the power. As the value of x/x. approaches 1, the first order approximation 
of PMD fails and equation 10 is no longer valid. 

Referring to Fig. 8, the linear polarization states represented by circle 570 are 
ao synonymous with equator 510 of tire Poincar. sphere 500 (see Fig. 4). Assuming that optrcal 
fiber 22 includes an x-horizontal linear PSP 575 and the y-vertical linear PSP 577, ,.e., two 
orihogonal PSPs, all otiter points on the circumference of circle 570 represent linear states at 
differem orientations. A SOP 578 represents one possible linear SOP of an optical srgn^ of 
optical fiber 22. SOP 578 contains components of both x and y polarizations, i.e., SOP 
,5 is a weighted linear combination of PSP 575 and PSP 577. Depending on the amount of t,me 
delay between PSP 575 and PSP 577, SOP 578 has a DOP tirat is less than or equal to untty. 
AS the angular distance, 2X. of SOP 578 to a PSP decreases, DOP increases. At tire cntical 
angular distance, a SOP 579 is equally distant ftom PSP 575 and PSP 577. i.e., SOP 
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579 is a 50/50 mixture of PSP 575 and PSP 577, resulting in a minimum DOP. The 
definitions of the Stokes parameters, such as S„ provide a relation between a and the radial 
distance IX such that Equation 10 takes the form. 



DOP(t,2A.) = 



1 + 



-2 




cos 



(2?i) 



(11) 



In general, the DOP is a fiinction of 2^, the angular distance between the SOP and any PSP 
on the Poincare sphere. IV is a function of the longitude distance, IX, and the latitude 
distance, 2w. Using the definitions of the Stokes parameters. Equation 10 is rewritten in the 
form. 



D0P(t,2?i') = 



1 + 



^-1 



"^0 V'^o 



1 + 



-1-2 



{2V) 



cos 



(12) 
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If the PSPs of an optical fiber are unknown, computer 120 runs algorithm 600 while 
polarization modulator 20 changes the SOP. For each SOP, algorithm 600 calculates both 
the location of the SOP on the Poincare sphere and the DOP. The algorithm sends these 
values into memory and repeats the cycle. Algorithm 600 stops collecting data points and 
fits e g , by linear-least-squares, the data in memory to find the maximum DOP, i.e., a DOP 
of unity corresponds to a SOP which represents one of the PSPs. Algorithm 600 collects 
enough data points so that the DOP as a fimction of SOP is well represented. Collecting too 
few data points leads to incorrect fitting results, whereas collecting too many data points is 
time consuming and allows enviromnental changes, i.e., temperature and stress on the fiber, 
to affect the location of the PSPs on the Poincare sphere. Once the CPU calculates the 
identity of the PSPs, algorithm 600 calculates a polarization transformation necessary to 
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transform the PSPs into linearly polarized x-horizontal and y-vertical PSPs aligned to the x- 
horizontal and y-vertical axis of delay module 170. 

Algorithm 600 uses Stokes parameters and Jones matrices representmg the PSP 
polarization states, retarder 140, and retarder 150 to calculate the settings of the retarders 
which transform the PSPs into linear horizontal and vertical states. Algorithm 600 begms 
with an arbitrary polarization state A 



15 



20 



(13) 



where 5^ + = U and a horizontal polarization state, i.e. a x-horizontal state, is represented 
1 0 by the Jones matrix, 



X - horizontal = 



(14) 



25 



The angle. A, between a major axis of the arbitrary state and the axis of a Cartesian 
coordinate system is 

A = 2tan-'[(2 5 8cos((p))/8'-e^)]. 0^) 
Algorithm 600 converts the known PSPs from Cartesian coordinates into Jones matrices, and 

determines A. 

Once A is known, Algorithm 600 calculates how to orient retarder 140 to convert the 
PSPs into linearly polarized states having an angle. A, between the x-axis of the Cartesian 
coordinate system and the axis of the linearly polarized states. Next, Algorithm 600 
calculates how to orient retarder 1 50 to rotate, i.e., by P, the linearly polarized states so that 
they coincide with x-horizontal and y-vertical polarization states. 

A complete transformation of an arbitrary state into a linearly x-horizontal 
polarization state using a quarter-waveplate as retarder 1 40 and a half-waveplate as retarder 
150 is 

[R(-p/2) X HW X R(P/2)1 X [R(-A) X QW X R(A)] (16) 



where 
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R{D) = 



cos(D) sin(D) 
- sin(D) cos(I>)_ 



(17) 



HW = 



i 0 
0 -i 



and 



(18) 



QW = 



0 



0 



(19) 



The computer sends electrical control signals, based on the transformation 
calculations, both to retarder 140 to transform the PSPs to linear PSPs and to retarder 1 50 to 
rotate the linear PSPs to the x-horizontal and y-horizontal axis of the delay module. 

The transformed optical signal propagates into delay module 170. Referring back to 
Fig 3, polarizing beam splitter 171 redirects x-horizontally polarized optical signal, i.e., PSP 
2 through a variable delay line and y-vertically polarized optical signal, i.e., PSP 1, through 
a fixed delay line. Polarizing beam splitter 1 72 recombines the x-horizontally and y- 
vertically polarized optical signals. Before the recombined signal exits delay controller 200, 
beam splitter 165 redirects a portion of the optical signal into polarimeter 210. Polanmeter 
210 is similar to polarimeter 110 described above and sends electrical signal proportional to 
the stokes parameters to control circuit 220. 

The control circuit, e.g., a microprocessor, calculates the DOP of the recombined 
signal AS shown in Equation 10, DOP is a function of both the time delay, x, and the SOP of 
the optical signal, i.e., the relative power, a, of each PSP. Control circuit 220 calculates 
DOP and determines a time average DOP as polarization modulator 20 modulates the SOP of 
the input optical signal and the relative power, a, of each PSP. Referring back to Fig. 7b, the 
time averaged DOP, i.e., as a changes between values of 0 and 1, is at maximum for zero 
delay between PSPl and PSP2, whereas the time averaged DOP decreases as the delay 
between the PSPs increases. Control circuit 220 sends electrical signals to delay assembly 
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174 to adjust the time delay between PSPl and PSP2 sueh that the trnte averaged DOP is 

maximized. 

Control circuit 220 averages the DOP for a time period that is sufficient to 
characterize the DOP for several SOPs. The shortest DOP averaging time period is set by the 
rate at which polarization controller 100 aligns the PSPs to the x- and y-axis of the delay 
module. At longer time periods, collecting too many data points is time consuming and 
allows environmental changes, i.e.. temperature and stress on the fib«, to affect dre locaUon 

of the PSPs on the Poincare sphere. 

in an alternative embodiment, referring to Fig, 9, a CPU 720 of a compensating 
apparatus 700 receives electrical signals ftom a polarimeter 730 and a control unit and runs 
an algorithm, as wiU be described below, to determine the location of the two PSPs of optical 
fiber 740 on the Poincare sphere, and the magnitude of PMD time delay between and the 

relative order, in time, of.be two PSPs. We Control unit includes a DOP sensor 770. such a. 
a polarimeter or a spectrum analyzer, and a CPU 780. DOP sensor 770 monitors a po«»n of 
the optical sign^ exiting a delay module 760 and sends electric signals to CPU 780. CPU 
780. in turn, calculates a value of DOP. DOPs,™. based on *e electric sign^s sent ftom 
DOP sensor 770 and sends an electrical signal proportional to DOPs,oN back to CPU 720^ 
CPU 720, in turn, uses this value of DOPs,™ to determine the relative order, in tmre, of the 
two orthogonal PSPs, e.g., which PSP leads and which PSP trails. 

During operation, CPU 720 determmes which PSP is leading and which >s tradmg by 
comparing two different DOPsioN values, one DOPs,™ value determined when the delay 
module retards one PSP and the other DOPs,™ value detemiined when the delay module 
retards the other PSP. A higher DOPs,oN value corresponds to less PMD time delay and a 
correctly retarded PSP. When the wrong PSP is delayed, i.e.. the trailing PSP .s further 
5 delayed in time relative to the leading PSP. the PMD time delay increases and the DOPsio. 
value decreases. Typically. CPU 720 records »d compares two values of DOPs-oN at a t^e 
faster than the rate at which PMD changes in the optical fiber. After determining whrch PSP 
is leading and which is trailing. CPU 720 continuously monitors the DOPsioN value 
calculated by CPU 780 to determine whether or not the leading and trailing PSPs have 
„ switched order. If the PSP order switches, the DOPs,o. value will decrease because Are 
delay module will further delay the "new" trailing PSP the^by generating more PMD t,me 
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delay be«vee„ .he two PSPs. When CPU 720 a decrease in the DOPs™ value, the 

CPU uses this information to generate electric signals which when sent to a polanzatton 
transformer 750 cause the transformer to rotate U>e PSPs relative to the delay module such 
that the delay module will retard the -new" trailing PSP. Alternatively, without rotatmg the 
polarization of the PSPs. CPU 720 can send an electric signal to the delay module caustng 
the module to change which PSP U retarded, such as by changmg the optical path delay. See 

CPU 720 runs an algorithm that records ele*ic signals ftom the polarimeter and the 
control unit and stores these values into memory. The algorithm processes d>e stored values 
from the polarimeter to determine where the PSPs of the optical fiber are located on «>e 
Poincar. sphere and the magnitude of the time separating flre two PSPs. The algortthm also 
processes the stored values from the cont^l unit, as discussed above, to determme the 
relative order of the PSPs. Once CPU 720 has detemtined the location of the PSPs and 
identified the leading and trailmg PSP, CPU 720 calculates a polarizafion transforma.,^ 

will orient the PSPs relative to a delay module 760 such that the module will retard the 
leading PSP relative to the trailing PSP. CPU 720. in turn, sends electrical signals to 
polarization transformer 750 to produce the calculated polarization fransformahon. CPU 720 
determines which electric signals will generate the desired transformation by using a lookup 
table containhrg input polarizations, values of applied electric signal — 
polarizations. Tlte lookup table can be preprogrammed by a system user mto CPU 720 pnor 
„ operation or determined via a calibration procedure described below After determtnmg 
the magnitude of the PMD time delay. CPU 720 sends electrical signals to delay module 760 
to set the overall magnitude of time delay between the two PSPs. 

CPU 270 processes each of the signals ftom the polarimeter to generate a fme 
averaged state otpolarization of a stokes vector. 7 . CPU 270 uses the time averaged state of 
polarization to extract an estimate of the location otthe PSPs and fte magmh.de otthen 
separation via a pertubative expansion of r . 

where 7 is the time averaged state of polarization of a stokes vector measuredby the 
30 polarimeter. 's ... and 7 .,„ a. two components otthe SOP. - *e signal bandwidth, 

-16- 



25 



Docket No.: 01997-296001 

>l2 J m 



0 ^ 

is ,he magnitude of the PMD time delay and can be expressed as |nf , and 0! W} 
„.s higher-order terms of the pertubative expansion. Equation 20, minus the htgher- 
order terms. 0(A»n. is equivalent to equation .0 and is generated by starting witir a two- 
mode, pure-polarization input field, written .(,). a(,) is transformed by an optical fiber mto 
5 an output field 6^?;, expressed as 

ft(0 = r(w)a(0=jf/(^)^""nw)«o (2^^ 
where/is the signal spectnun and Jis a complex 2-din.ensional vector. T^. time averaged 
state of polarization is represented by 

r = |dtb'(t)ob(t) P2) 

,0 where S is a vector of Pauli matrices. Using equations 21 and 22, the time averaged state of 

polarization can be expressed as 

PMD is chamoterized by a (Stokes-space) PMD vector, 5 , which can be defined by the 

relation 

— = Q(0))xs(c0) (24) 

^5 dco 

By using equation 24, and by defining a center of the signal spectrum coo as 

J In 

and a nominal SOP as 

.0 a pcrturbative expansion for time averaged state of polarization can be expressed as 

Q.x;..aK(5.1.)]+0{A.,'} (25) 



r - So- 



■ • 2 L 

where ti,e signal bandwidth is characterized by Aa,^. ' ^ H'lCw-w.)' /2;r . 
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The Q X a X term can be described physically as a reduction of components orthogonal to Q 
due to first-order PMD and the Q^x so term can be described physically as a rotation due to 
second-order PMD. Rearranging equation 25 yields 

- o ^ Ar,^2 s I +-Aco^ Q xs„+o{aco' (26) 

Using equation 26, the physical effects of PMD of the optical fiber on an SOP can be 
described physically as generating a rotation and a rescaling of the SOP: 

r =[Rotate {Clj] [scale{Q.)\s,. (^7) 
In terms of first-order PMD, a rescaling means that DOP is a function of the angle 
between and Q . Although discussed in terms of first- and second- order PMD, the 
perturbative expansion for time averaged state of polarization can be used to develop 
analytical models for higher-order PMD which can be used to extract higher-order PMD 
parameters from DOP measurements. For instance, a CPU in a PMD compensating 
apparatus capable of compensating higher-order PMD can use the higher-order PMD 
parameters derived from equation 20 to generate control signals to reduce higher-order PMD 

15 in the optical signal. 

For simplicity, we limit our discussion to first-order PMD calculations. Given a 
number, N, of time-averaged SOP measurements r, n , resulting from unknown and 
varying input polarization input fields ao, the samples, n ideally obey the relation 

(28) 

20 for first-order PMD, where Z is an unknown vector to be estimated by the CPU. In Stokes 
space, forms the surface of an ellipsoid within the Poincare sphere whose major axis is 

aligned with the PSP directions +a and - a. Vector 2 satisfies the relationship (1-| a = 1 - 
Aco^„.se x^n. which relates the magnitude of 2 to the PMD time delay between the PSPs. 
The CPU can estimate the vector a by mathematical methods such as eigenvalue 
25 calculations or nonlinear curve-fitting. In nonlinear curve-fitting, the CPU records a series of 



DOP data, r*, and coverts them into magnitudes, dk □ 



rk 



and directions = r/ 



. The 
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CPU performs a nonlinear fit by minimizing the distance between the measured DOP (d,) 
andtheexpectedDOP(f(f., a)). For example, the least-squares estimate of a is 

5^,^,, =argminF,(fe) (29) 



15 



b 

where 



F,ib) = ^[d,-Mp)f ^^^^ 

k 

(31) 



'^'^'^ = K-(b-r.)b| 
where the function f(f„ b ) gives the theoretical degree of polarization (DOP), or 

magnitude of r . 

In eigenvalue calculations, the CPU can estimate a by using matrix method 
10 techniques to maximize the second moment of the points on an unknown projection. 

Typically, CPU rescales the measured DOPs, r „ before performing moments calculations 
to artificially sharpen the contrast in DOP. CPU rescales the time-averaged SOP 
measurements via the relationship 

(32) 

q, =S(|fJ)f, 

. . //T^rMi T\r\\) ^/1 nop ^ CPU uses the scaled vectors, qic, to 
where S(lrk|) is equal to ((DOPk-DOPmin)/l-UUi'm.nj. u ubcb 

calculate an estimate of the PSP axis 

=argminF2(6) (33) 

1^1=1 

where 

The funcln generally reaehes a maximum near the direction of fte PSPs, because 
PMD tends .0 pull the SOP towards this axis. CPU can also use the unsealed time-averaged 
DOP measurements, T , to estimate I , but large values of DOP„. decrease the contrast 
between the measurements which, in turn, can decrease the precision otthe method because 



20 
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.andom „o„™fornu«es in the input po,a,iza«o„s can domimte ^ Motion F. and reduce 
d,e accuracy of detennming I Equation 34 can be r^ged to provide 



10 



15 



^2 = ^<l/j 



V k 



b^b^Mb 



(35) 



which reduces to the very simple form: 

equation 36 canl. solved using Lagrange multipliers. Tlre estimated , is the eig— 
corresponding to the largest eigenvalue of the 3x3 matrix M. Hre magmtude of the PMD 
tinte delay is equivalent to flre magnitude of the vector I . The CPU determines the 
.nagnitude of I by analyzing the statistics, such as mean and minimum, of the DOP data. 

Random flucU>a,ions in the density of sampled SOP's can effect the precision at 
which CPU determines the estimated 1 . To increase precision, the CPU can record more 
samples. Recording more samples reduces *e overall speed a. which d,e compensaUng 
apparatus can reduce PMD in the optical signal. Alternatively, the CPU can increase 
precision by correcting random flucU^tions in the density of a sm^ler subset of record«l 
samples. CPU reduces the random fluctuation by using a moment function to remove DOP. 

V k J 

Tlre matrix method also has an advantage that much of the calculation can be 
performed in real-time as the CPU records data ftom dre polarimeter. Each data potnt can be 
used to increment a matrix. M. with very few floating operations by continuously updatmg M 
.0 to include new data and throw out old data. The CPU can update Musing ti,e ..lationshrp: 

menusingthis scheme, the CPU canmaintainanupdatedMat all times so that 
calculations, such as an eigenvalue decomposition to estimate the PSPs. can keep pace w.th 
data acquisition. 



-20- 



kDocketNo.: 01997-296001 



The vec». J , climated by CPU has degrees of freedom, one degree for a 
magni^deof |a| and «,o degrees for dre direcdon of t . T.e relationship be«,ee„ U,e 
magnitude of vector J and the PMD vector. 5, which also has dtxee degrees of fteedont, is 

n = ±(2|a|'/«.Ur^- 

where aW is 4^ "-^^"^ *^ "P""'' """"T 
:!::„ticalsign-band«idthdtat«asn,easuredprior.„operationhecansec^es.n 

PMD generally do no. affect bandwidth. AltemaUvely the compensating apparan. can 
ZndLspelmanalyzerwhichmonitorsthebandwidthoftheopticalst^andsends 

elect^nicsignalsproportionaltofltatbandwidthtotheCPU. 

When calculating an estimated PMD vector with equatron 39, the CPU rs unable o 
determine dtc correct Sign because the mathematical methods cannot be used to c«^^^^ 
direction either.aor-a,fromtheDOPdat.Asdiscussed above wrthres^ttoF^J^ 

compensating apparahrs 700 determines the correct sign by measurmg dre DOPs,on after dte 

Idulc in other embodiments, referring to Pig. 10, the CPU of a compensattng 

appLtus SOO determines the sign without monitoring a DOPs» value. CompensaUn 
:aratusS00incMesa,ixedPMDmodule810,suchasabire«ngen.crystalordeUy 

module shown in figure 3, which adds a fixed, known PMD time delay, n^^. to the 
variable, u,^ownPMDtimedelay,a«.cw».producedwi*intbeopdcal fiber. By ad m 

, .he fixed. ..own delay to the optical signal, CPU can .dentin the leading and trarlmg PS^ by 
Irmh.ngthec„rrectstgni„e<,uat.on39. CPU can determine the correct s.gn provtded 
,hataK.ow«isgreaterthana.K.ow..Forexamp.e.givcnarelati„nsh,psu^^^^^^^ 
(, 00.x)^ unambiguous values of x can be calculated for any measured value of z provrded 
drat X is less than 100. The exact value of axscw. depends upon d>e magnitude of 
a and can be estimated by a system operator prior to operation or can be determtned 

" ::::;:u iipi..hecpuc.^ ^ 

:r:yasc.ingfac.„rtoproduceavalueof.efixed,.owntimede,.d.^^^^^^^ 
relationship «kho« > "-tc The CPU, in mm, sends an eleComc srgnal to fixed PMD 
„>„dule 810 to adjust the fixed delay to produce dre necessary delay. 
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Referrtog to Fig. 1 1 , a system 900 for calibratmg a polarization transformer 940 
includes a transmiUer 905, a polarization modulator 9.0, a polarimeter 920, a CPU 930, and 
a calibrating module 945. As de^ribed above, during operation, CPU 930 sends conttol 
signals to polarization transformer 940 to transform an optical signal such that a delay 
5 module (no, shov™) reduces PMD by retarding the leading PSP Tl.e objective of system 900 
is ,0 determine, for a given set of ele^c control signals, what polarization state mput mto 
polarization transformer 940 will be tramformed to different polarization states defined by 
the calibrating module. Calibrating module 945 includes a polarizing beam sphtter (PBS) 
950 having an orientation p , and two photodetecti.rs 960 and 970. As a transformed SOP, 
,0 s from polarization tiansformer 940 p^pagates into calibrating module 945, PBS 950 
redirects different portions of s onto detectors 960 and 970 based on the polarization of s . 
Each detector produces separate currents tt.at are proportional to different portions of the 
ti-ansformed SOP The currents are given by the expressions 

(40) 

ji = 14(l+ s- p) 

,5 and ^^j^ 

i2='/.(l-S- p) 

Where /, and are normalized such ti,at (,-, . k) - 1 . An an^og device 946 ~ , and . 
and sends an electric signal proportional to tire difference of tire two currents, to CPU 930. 
U is the projection s ■ P . Alternatively, CPU 930 can receive /, and h and calculate ti,e 

ao difference i. As polarization modulator 910 changes the SOP of the optical signal 
transmitted by transmitter 905, CPU 930 determines the SOP of tire optical signal by 
analyzing the electiic signals from polarimeter 920. and records and ti,e magnitude of ti>e 
conttol signals being sem to polarization transformer 940 for that SOP By collecting a se, of 
SOP measurements, s, and storing them into a matrix S. CPU can consttuct an 

25 overdetermined maUix given by 

(42) 

Sp= id 

where each s;isa3-dimensio„alrowofS. p is a S-dimensional vector, and i, is tire vector 
of difference currents. CPU can estimate p by computing the pseudo-inverse of S. To 
calibrate polarization tiansformer 940, CPU compute an estimated p tor a plurality of 
electric control signals and SOPs. CPU stores each estimated p , SOP, and tite values of tire 
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electric control signals in a lookup table. The lookup table provides specific electric control 
signals that will align the SOP with the axis system of the PBS. When compensating PMD, 
the calibrating module is replaced with the delay module such that the relative orientation 
between the module and PBS is known. After determining the location of the PSPs, the 
magnitude of the time delay, and which PSP is leading, the CPU uses the lookup table and 
the known relative orientation between the module and the PBS to determine which electnc 
control signals are needed to transform the PSPs such that the delay module retards the 
leading PSP In other embodiments, referring to Fig. 12, a calibrating module 945 and a 
delay module, such as described above with respect to figure 3, can be integrated into a 
single unit 95 1 so that the CPU can create a lookup table specifically for the polarization 
transformer and the delay module being used in a compensating apparatus. 

In other embodiments, the polarization transformer can include polarization 
transforming devices, e.g., electrooptic, acoustooptic, or stress induce bifringence, which can 
transform the PSPs to linearly polarized PSPs aUgned with the x- and y-axes of the delay 
controller. The delay module also can include any optical device capable of delaying of one 
polarization state relative to another polarization state. Although shown aboveas separate 
units the polarization transformer and delay module can be integrated into a single unit that 
functions to transform the optical signal and delay the leading PSP to reduce PMD in the 
optical signal. An integrated polarization transformer and delay module is descnbed, for 
example, by C. Glingeneret. al. in Optical Fiber Communications Conference PD29, 1999. 
A measure of performance of the PMD compensating apparatuses can be determined by 
monitoring the jitter in the rising edge of the PMD compensated optical signal. A histogram 
of the rising edge includes a variance that is proportional to residual PMD. 

It should be understood that the foregoing description is intended to illustrate and not 
limit the scope of the invention, which is defined by the following claims. Other aspects, 
advantages, and modifications are within the scope of the following claims. 
What is claimed is: 



-23- 



